Quantum Joint Distributions * 



Todd A. Oliynyk t 
School of Mathematics and Statistics 
University of Canberra 
ACT 2601, Austraha 

February 1, 2008 



Abstract 

In this paper we provide a method for constructing joint distributions for an 
arbitrary set of observables on finite dimensional Hilbert spaces irrespective of 
whether the observables commute or not. These distributions have a number of 
desirable properties: they agree with the standard quantum mechanical ones if the 
observables commute, they also depend continuously on the observables, and under 
unitary transformations they behave in a reasonable manner. 

1 Introduction 

As discussed by Ballentine the assumption that a particle has a definite position 
and momentum is not incompatible with the uncertainty principle. More generally, the 
assumption that a quantum state has definite values for certain observables even if the 
observables do not commute is compatible with the statistical interpretation of quantum 
mechanics. Therefore, as indicated by Ballentine 0], there is, in principle, no problem 
with admitting joint distributions for any set of observables commutative or otherwise. 

As the standard quantum formalism does not provide a recipie for constructing joint 
distributions associated with non-commuting observables, the formalism must be ex- 
tended. The problem can be thought of as the problem of constructing joint distribu- 
tions with given marginals. In [5] (see also [HIS]), Cohen and Zaparovanny provided a 
method for constructing joint distributions with given marginals. It was later shown 
in jlUlll3) that the Cohen-Zaparovanny construction actually generated all possible joint 
distribtuions. Consequently, it is possible to construct joint distributions for any set of 
observables, commuting or not. The difficulty is that the Cohen- Zparovanny construction 
does not single out a particlular distribution. 

In this paper we provide a method for constructing joint distributions of an arbitrary 
set of observables on finite dimensional Hilbert spaces irrespective of whether the ob- 
servables commute or not. We restrict ourselves to observables with distinct eigenvalues. 
This is not entirely necessary because if there is some degeneracy then the construction 
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still goes through but does not result in a unique distribution. Instead, a number of 
distributions are possible and we do not have a criterion for selecting out a particular 
one. 

The method we employ is technically simple. We first represent the usual transition 
probabilities as the volume of certain regions in projective Hilbert space. The joint 
distributions are then constructed naturally via the volume of the intersected regions. 
This is the strength of our method as compared to 013111 because once we choose to 
represent the transition probablities as volumes then then the joint distributions arise 
naturally. For all the technical simplicity of our method, the drawback is that the joint 
distributions are difficult compute as they require the evaluation of complicated integrals 
on projective Hilbert space. However, we do show that these distributions have a number 
of desirable properties. For example, they agree with the standard ones if the observables 
commute. They also depend continuously on the observables. Finally, under unitary 
transformations they behave in a reasonable manner. 

We also note that the volume representation of transition probablities introduced here 
can be used to simplify the constructions used in the hidden measurment approach to 
quantum probabilites |2I3I6I5| . Moreover, the present approach to transition probabilites 
can also be used to easily construct new hidden measurement systems. We will report 
on this work elsewhere. 

2 Projective Hilbert Space 

In this section we review some basic results about projective Hilbert space. We use the 
book jJT] as our standard reference. Let (H, (•!•)) be a complex Hilbert space where the 
inner product (-j-) is taken to be linear in the second variable. Define 

:= H \ {0} and := C \ {0}. 

On Ti.^ we can define an equivalence relation ~ by 

ip ^ (j) a and only if there exists a A G such that = X(j). 

Letting [ip] denote the equivalence class for -0 e 7i^, we have 

[^] ■.= {X^/j\XeC''} . 

Projective Hilbert FH space is defined as 

FH := W / ^= {[ip]]^; £ W } . 

It is well known that P7i carries a Hilbert manifold structure for which the canonical 
projection 

is a C°° submersion. As a consequence for any q e FH and Vq G TgFH there exists a 
ip £ and (p £ H such that 

q — \>p\ and Vq — T^tt • (p . 
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This can be used to define a complex structure J and a strongly non-degenerate symplectic 
form u) on P?i via the formulas 

J(T^7r •(/)):= T^, • (2.1) 

and 

u;[^](T^7r • 01, T^TT • ^2) := 2;i|lV'|r^Im ((0i|02>||Vf " (<^i|V')(V'l02)) (2.2) 
for every xp G and (/>, ^i, 02 € H. It should be noted that 

givjw) — uj{v,^w) (2.3) 

defines a Riemannian metric on P7i and hence establishes that VTi is a Kahler manifold. 

Given a function / e C°°(PW), the fact that the symplectic form is strongly non- 
degenerate implies that the following equation 

ixfUJ = df (2.4) 

uniquely defines the vector field Xf. The Poisson bracket {•, •} is then defined via 

{f,g}:=iv{Xf,Xf) y f,geC'^{rH) . (2.5) 

Let C{H) denote the set of bounded linear operators on H. Then the unitary group 
U(W) is defined by 

U(W) := {U e £(W) I {Uil>\U^) = VV-, e W } . 

Its Lie algebra u{Ti.) is the set of skew-adjoint operators, i.e. 

u{n) -.^{Ae c{n)\A^ = -A}. 

Here we are using f to denote the adjoint of an operator. The following map 

p : V{H) xFH-^PH : ([/, [tp]) [Uip] 

defines an action of U(H) on PTi by symplectomorphism (i.e. p^w = u) for all U £ U(H)). 
There exists an equivariant momentum mapping J : FH — > u{H)* for this action defined 

by 

(J([V]),^) :=in%^ V^eW, Aeu(W), (2.6) 



where (•,•) denotes the canonical pairing between u(W)* and u(W). Letting C°°(PW) 
denote the set of smooth functions on FH, the momentum map can be viewed as a map 
J : u{H) — > C°°(PH) by defining 

J(^)(a;) := (J(a;), A) V a; e PH . (2.7) 

Because of the equivariance of the momentum map J, we have the useful property 

{J{A),J{B)} = J{[A,B]) \fA,BGu{n). (2.8) 

Let Sa denote the set of bounded self-adjoint operators on H. For each operator 
i? e Sa we define a smooth function {H) on FH via 

{H):=3{-iH). (2.9) 
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Thus {H) is just the usual expectation of the observable H , i.e. 

W(l*l) = ^- ("0, 

With this notation H2.8|) can be written in the more familiar form 

{(A),(i3)} = (l[S,A]). (2.11) 

3 Action angle coordinates on fH 

For the remainder of this article, we will assume that diniTi = iV < oo. Let {-01, • ■ ■ , V'Ar} 
be an orthonormal basis for 7i. Define the projection operators 

P^,^Hk){M- (3.1) 
We can use the momentum map to define smooth functions p^^. on VH by 

P^.= {P^,). (3.2) 

Using H2.10|l we get 

which is the transition probability from the state ip to ipk- As the operators P^^, commute, 
formula H2.11|l shows that the functions {p^^, . . . are in involution, i.e. 

{V„V'fc} = Vj,A; = l,2,...iV. (3.4) 
It follows from I = -Pif>k that 

N 

X^PV-fc^l' (3-5) 

which shows that at most {N — 1) of the functions p^^. can be independent. It is not 
hard to show that the set {p^2^ ■ ■ ■ jPV'n} independent. That is the set of points in 
IPT^ for which the covectors {dp^^, . . . are linearly dependent has measure zero. 

Consequently, we can use these functions to construct action angle coordinates of PH. 
following the standard recipe, see for details. This results in the following coordinate 
chart 

T : T^-'^ X S — >Fn 

N N 

{e, I) = {{02, . . . , 9n), ih, ■■■■Jn))^ [(1 - J2 + E ^"''^ VIjM (^-^^ 

fe=2 j=2 

where T^~^ is the {N — 1) torus and 

N 

5:= {(/2,...,/iv) eM^-MO</, ^/, <l}. (3.7) 

i=2 
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In this chart, the symplectic form lo is given be 

N 



uj = hJ2dOk^dIk. (3.8) 

i=2 



We also note that the functions p^^ have the coordinate representations 

p^^i9J)^Ij j = 2,3,...,iV. 
We can define a volume form fi on PTi. by 

/^-(^j i^A^. (3.9) 

N-l times 

Locally this given by 

(N - 1)! 

M = (27r)(^'-i) A . . . A 6lAr A (i/2 A . . . A d/AT . (3.10) 

Then because the chart ()3.6|) covers all of P7i except for a set of measure zero, the volume 
of PH is given by 

Vol(PH) = -^^^^^^ ...^ de2...d9N j^dh...dlN. 

A straightforward calculation shows that 

/ dh...dlN = ll{N 
Js 

and hence Vol(P7i:) = 1. 
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4 Volume representation of transition probabilities 



Suppose ip^ (j) Cz Ti^ . Then the transition probability from the state ^ to 0, or vice versa, 
is given by 

urn? 

Because this formula is invariant under scaling of or by non-zero complex numbers, 
it passes to a well defined function on FH x FH given by 



m'/'):=^^:?f^- (4.1) 



It was shown in that if we let d{x, y) denote the geodesic distance between points 
X, y £ FTi then the distance d{x, y) is related to the transition probability T{x, y) via the 
formula 

T(x,y) = cos2(^). (4.3) 
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Figure 1: Volume representation of transition probabilities 

This shows that there exists a representation of the transition probability in terms of the 
geodesic distance. The question now is, are there other representations for the transition 
probability in terms of geometrical objects on VH1 We will show that the transition 
probability, at least for finite dimensional Hilbert spaces, can be related the volume 
of certain regions in FH. To motivate this, we will first look at FH where 7i is a 2 
dimensional Hilbert space. Then FH = S*^ where S*^ is the ordinary two sphere in M.^^. 
Suppose {tpi-, 4'2} is an orthonormal basis for Ti. and ip € Ti.^ is an arbitrary state vector. 
Since {^pl^^p2} is orthonormal, we can choose them to be the north and south poles of 
S'^ as in Figure ^ The symplectic form on is 

^sin6' , , 

u; = dm A dO 

2 

while the volume form /i is given by 

The normalization on the volume form is chosen so that fi = 1. 

Referring to Figure ^ let fl be the shaded region between the points [-01] and 
Then a straightforward calculation shows that 

r([V'2],[V'])-Voi(r!) /m. 

Of course if we let 7 be the geodesic between [^p] and [i/;2] represented by the dashed line 
in Figure n then we also have 

T([V'2], ^ cos2((2n)-\geodesic lenth of 7)) . 
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Letting fl'^ denote the complement of ft we also have 

r([v.i],[^]) = Voi(r!^) := [ ^i. 

It is interesting to note that the conservation of probability 

i = r([^i],[^]) + T([^2],M) 

has the simple geometric representation 

¥n = nun'' . 

Now, consider the self-adjoint operator 

A = ai|Vi)(i/'i|+a2|V2)(i/'2| 

where the eigenvalues ai and 02 are distinct. Then upon measurement, the probability 
of a measurement of the observable A yielding a value aj given that it is in the state [ip] 
is 

Prob(^ = a, I [^]) - Vol(%) J - 1, 2 , 

where we let 

^2 := and fli := . 
Now suppose {(/>i, (f>2} is another orthonormal basis and B is the operator, 

B ^bi\^l)((j)l\+b2\(b2){h\- 

Then as above, we find two sets Si and E2 such that 

PH = Si U S2 

and 

Prob(B = bj I [^A]) = Vol(Sj ) j = 1, 2 . 

So far we have not done anything new, except express the usual transition probabilities 
in terms of the volume of certain regions in P7Y. However, we observe that there is a 
natural way to construct a joint distribution for the observables A and B which is given 

by 

Prob(A = aj,B = bk \ [V']) := Vol(f7j n Efc) j, /c = 1, 2 . 

It is important to realize that this formula is valid even if A and B do not commute. It 
is also easy to see that it gives the correct joint distribution if A and B do commute. See 
figure 121 for an illustration of the region ^2 r\Y,2, the volume of which would yield the 
transition probability Prob(y4 = 02, S = 62 I [V'])- 
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Figure 2: Volume representation of joint probabilities 

5 Joint distributions 

To generalize the construction in the previous section to arbitrary but finite dimensions 
we must first find a method for generalizing the decomposition FH = f^i ur22- So for the 
moment, let us still assume that dimW = 2 and that {ipi,tp2} is an orthonormal basis. 
Now, a short calculation shows that 

X GClk if and only if p^. {x)p^^ {y) > p^^ {x)p^^ {y) for j = 1, 2 , 

where 

y = [ip]. 

This motivates us to make the following definition. Let 7i be an iV dimensional Hilbert 
space. Suppose /3 = {tpi, ■ ■ ■ , ^^n} is an orthonormal basis for H.. Then define a region 
fl{y,f3,tpj) of FH that depends on a point y e PH, the basis (3, and a particular basis 
vector ipj by 

^{y, P,ipk) ■■= {x gVH\ p^^ {x)p^^ (y) > pv-fc {x)prpj (y) for j = 1, 2, . . . , iV } . 

It is useful to introduce an alternate characterization for /3, t/i^) which although 
it looks more complicated is actually easier to work with. To start, consider the following 
vectors in R^~^ 

jth 

a-0 and ^,+1 := (0,..., 1,...,0) j = 1, 2, . . . , TV - 1 . 

Let 

N 

V-=Y.p^M^^ (5-1) 
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and define 

N N 

{ J2 lj^i+hi^^^'^\lk>0 and ^/j-l}- (5.2) 



Proposition 5.1. Suppose (3 = {V'li ■ • ■ iV'-Zv} *s an orthonormal basis for TL andy G P7i. 

N 

Proof. Assume that p^/,^ (y) 7^ 0. The case p^^. (y) = wUl be left to the reader. Then 
using 1)5. we can write ^^^=1 Pipj^^)^] 

S ^ y'"' ~ ^f'- V + ■ 

From H3.5|l it is easy to see that 

N 

(5.3) 



2^ P^, a; ^Pi,, (x) + ^ = 1 



These two results along with the fact that p^. > show by the definition of S{y, f3, ipk) 
that 

N-l 

^PM'oi^)^] ^ S{y,(3,^k) P'4,,{x)p.4,^{y)>p^^{y)p^^{x) j = 1, 2, 3, . . . , iV . 

□ 

The next two propositions show that the sets ^{y, /3, ipj) have the required properties 
to be considered a generalization of the sets fii and ^2 from the previous section. 

Proposition 5.2. Suppose 13 = {V'l, ■ • ■ , V'w} orthonormal basis forH and y e P7i. 
Then 

N 

Fn=[j n{y,p,i;,) 

fc=l 

and 

Yo\{n{y,(3,i;,)nn{y,P,^k))^0 forj^k. 
Proof. Let S denote the closure of S defined by H3.7|l. i.e. 

AT 

S = |(a2, . . . , aw) G M^"^ I % > and ^ < 1 | 
and define a map J : ^ R^^^ by 

N 

i=2 
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Since p^- > and X^jLiPVj 1' have that J(P'W) C S. To sec that this inclusion is 
actually an equaHty, consider any state vector 



N 
J = l 



where at least one of the coefficients Cj is non-zero. Then G 7i ^ and 

1 ^ 

It follows directly from this formula that i{WH) = S. Also, it is not hard to verify that 

N 

i=i 

The above two results and proposition 15 . II then imply that YH = UjLi ^{v-, Pi V'i)- 

From proposition lS.ll and the definition of J, we have that /3, ijjj) = J {S{y^ 13, ipj))- 
Consequently 

But for j 7^ fc, the set S{y,(3,4'j) ^ S{y,f3,4'k) lies inside an iV — 2 dimensional subset 
of R^^^ and hence J iS{y, (3, ipj) n S{y, (3, ij^k)) must have measure zero. Therefore for 
j T^k the formula Vo\{Vt{y, /3, i),) n 17(2/, V-fc)) = follows. □ 

Proposition 5.3. Suppose {ipi, . . . ,i^N} is an orthonormal basis for Ti and y G P7Y. 

Then 

Yo\{n{y, P, i^,)) = p^^ (y) J = 1, 2, . . . , TV . 

Proof. It is enough to prove it for j = N. From proposition 15. II and equations 13. 8|) - 
H3.10|l it is clear that 



Vo\{n{y,(3,M)= / / ^de^-.-dON [ dh...dlN 

= (iV-1)! / dl2...dlN 
But from (|5.1|) and (|5.2|) it is easy to verify that 

[ dh... diN = ttt^-tttpvjv (y) > 

which completes the proof. □ 

We are now ready to construct the joint distribution. We only handle the case of two 
observables with the generalization to more observables being obvious. 
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Theorem 5.4. Suppose A and B are two self adjoint operators with distinct eigenvalues. 
Let A — X^jLi '^j IV"/) ("0/1 '^nd B = X^jLi ^j\i^f){'4'f\ be the unique spectral representa- 
tions of A and B. Then for each y g WH, 

Prob(A - a„B^ bk\y) := Vo\{n{y , f]"" , ^f) n n{y,P^,^pf)) 

is a well defined probability distribution. Moreover, the marginal probability distributions 

N 

Prob(A = aj\y) := ^ Prob(A = aj,B = bk\y) 
fc=i 

and 

N 

Pmh{B = bk\y) := ^Prob(A = a„ S = bk\y) 

satisfy 

Prob(^ = aj\y) = p^,A{y) and Prob(B = bk\y) = P^b (y) , 
and hence agree with the usual quantum probability distributions. 
Proof. From proposition 15.21 we get 

N N 

j=i k=i 

N 

Also from proposition 15. 21 we get that 



Vol 



r Vol(O(y,/3^,V^/)nr!(y,/3-4,v.f4))=0 ii^^l 
- \ VolMj/,/3^>f)nf^(2/,/3^,0)=.0 iffc^m ' ^^'^^ 

Combining these two results shows that 

N 

l = Vol(PK)= ^ Vol(r!(2/,/3-4,V;/)nl7(y,/3^>f)), 

and hence 1 = "Y^^k^i Prob(v4 — aj,B — bk\y) ■ This along with the fact that Prob(A = 
aj,B^ bk\y) = Yol{n{y,P^,i,f )nniy,P^,i;J^)) > shows that Prob(A = 0^,5 = bk\y) 
is a well defined probability distribution. 
Now, 

JV 



= |Jf2(y,/3^,V/)nf7(y,/3^>f). 



fc=i 
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This result along with (|5.4|) shows that 



JV 

Vol(r!(y, /3^, ^/)) = J2 Vol(^^(y, /3^, ^f) n n{y, /3^, )) . 
fc=i 

Using this and proposition we find that 

AT 

^ Pi'ob(Aj = oj, Bfc ^ bk\y) , 

k=l 

as required. The same arguments can be used to show that 

Pi,^iy) = Y.f=i Prob(Aj = aj,Bk = bk\y). □ 

It is important to reahze that the above proposition is also valid for operators without 
distinct eigenvalues once a basis of orthonormal eigenvectors is fixed. The problem is that 
unhke the case where the eigenvalues are distinct there is no unique choice of basis and 
it is not hard to see that the joint distribution will depend on the basis that is chosen. 
Therefore, two different bases will, in general, lead to two different joint distributions. 
What is needed is a reasonable criteria that will select a particular basis. 

The next proposition shows that the joint distributions transform correctly under 
unitary transformations. 

Proposition 5.5. Suppose U e U{U), y G PH, and A = Ej^i and B = 

X]j=i IV"^) (V'^ I are two self-adjoint operators with distinct eigenvalues. Let 

N 

A' := UAU-^ = ^ a,\^f'){^f' \ i,f U^f , 

JV 

B' := UBU-' = J2 )(V'f I U^jf . 



Then 

Prob(v4' - a,,B' = bk\y) - Prob(A ^aj,B^ bk\puiy)) ■ 

Proof. To start, suppose a := {V'li ■ ■ • j V'Jv} is an orthonormal basis. From the definition 
of the functions it is easy to see that 

pui,, (y) = PV'j (puiy)) i = 1, 2, . . . , iv . 

It then follows directly from the definition of the set il.{y, a, ipj) that 

^{pu{y),a,ipj) = pu{^{y,a',tpj)) 
where a' :— := Uipi, . . . ,ip'j^ :— C/i/'jv}- Therefore, 

n{pu{y),f3^,^pf) n n{pu{y),f3'',i^i) = Pu {n{y, f3^' ,i>f) n n{y, V^f )) (5.5) 
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for all j,k — 1,2, N. Since p^jcu = lu it follows that the volume form /i is invariant 
under the action of U , i.e. p^fl = /i. So H5.5|l and the change of variables theorem implies 
that 

Yo\{n{pu{y),f3^, i^f ) n (y), )) = Yo\{^l{y, n^' ,i,f) n ^{y, v^f ')) 

for all j, A: = 1, 2, . . . , iV, which completes the proof. □ 

6 Continuity of the joint distributions 

Let dn denote the Hausdorff distance, i.e. if S C P7Y then 

dniA, B) ~ inf{ e\Ac B{B, e), B C B{A, e)} 

where for any set C C PH we define 

B(C,e) := {x erH\d{x,C) < e} . 

Proposition 6.1. Let [3 = {ipi, ■ ■ ■ , V'Jv} 7 ~ {^i, . . . , (p^} be two orthonormal bases 
for Ti. and suppose x £ P7Y. Then for any e > there exists a S > such that 

da(17(a;,/3,V,),r!(x,7,0j))<e j = l,2,...,N 

whenever — ipj\\ < S for j = 1,2, . . . , N . 

Proof. For j ^ k, 

I i^Mk) I = I (V^, - <PMk) I < 11^, ^<i>MM\<5i = 5. (6.1) 

Without loss of generality assume 5 <1. Then ||'0j ~ '/'ill < ^ implies that 

|2-(^,|0,)-(<^,|^,)| <<5, (6.2) 

and 

\m<i^3) - ('/'.iV'.)i = KV', - <i^M3) ~ {'PM, ~ '^,)i < A'pj - n\ u-j\\ < 2^ ■ (6-3) 

Let ip E H satisfy \\tlj\\ = 1 and x — We will assume that there exists a > 
such that 

|(V'|V,)I, l(V'l'/',)l J = 1,2,... TV. 

The case where |('0|'0j)l = or |(V'|0j)| = for some j can be taken care of by a slight 
modification of the arguments below and will be left to the reader. 

Fix e > and k € {1, 2, . . . , TV}. Suppose (j) e Ti satisfies ||(/>|| = 1 and [0] S 
n(x,7, ^ife). Then by the definition of n{x,j,4>k), ( [</>] )P0fc ( ["0] ) > P<l>ki['t'])P<PjibP]) for 
j = 1, 2, . . . A^, or equivalently 

\{ci^m^\iM4')\'>\{Am^mm' j^i,2,...n. (6.4) 

Define 

N 

vI/:=^(0,|0)A,V, 

3 = 1 
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where 



A, 



I I V') I 



Then for j 1, 2, . . . TV, 

and hence \^] G Q{x, P,^pk)- 

Note that Xj = KV^j - + 1| and so 

,\{t^^Ml<x <i , 



by 



But 



and hence 



mm 



S S 
1--<A,<1 + - j = l,2,...,iV. 



Now, ||vl/||2 ^^.^^ 7V|(0,|0)|2A2 and 1 - ||0|P = K'/'jl'^)!' imphes that 



1-J)<|l*f<(l + ^ 



(6.5) 



(6.6) 



by H6.5() . Also we have 



TV TV 
N N 

AT W N 

= |Ei('^i'^^)i'-Ei('^i'^^)i'(i-^^(^^-i'^^-))| - E \m,)^jm^i){m\ 

TV TV TV 

>|Ei('^i'^^)i1-Ei('^i'^^-)i'i(i-^^-(^^i'^^-))i- E i(0i0.>A,(^,#/)(</';i0)i. 



But note that 



TV 



(6.7) 



(6.8) 
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and 



iA,(V',i0,) - 11- i(A, - 1 + - ii < |A, - iii(^,#,>i + mm - ii 

< |A, - 1111^,11110,11 + Ik^jMj) - 1 + Wj) - i| + Imcl^,) - (0,l^,>l 

<- + Us-(^+'^S (6.9) 



where in deriving the last inequaUty we used H6.2() . (|6.3|l . and (|6.5() . Using H6.1|l . (|6.5() . 
H6.8|l . and (|6.9|) . we can write (|6.7|) as 

\{m\ > 1 - TV (^i + <5 - (TV - 1)2 J + 1^ 5 . (6.10) 

Therefore (jHinil and ^^U^ imply that 

|(vl/|0)p (l-jV(l/^ + 3/2)J- (TV -1)2(^^ + 1)^)2 
- ll^Pll^ilP - ((5/;. + 1)2 

or by 1121 and 



1 > cos2 



f dm,[m > (1 - iv(i/^ + 3/2)^ - (TV - ins/iy + 1)^)^ 
I %/2?^ y- (5/^+1)2 



This shows that by choosing S > small enough we get d([^],[0]) < e. Since [^] S 
J7(a;, /?, V'/c) and [0] € r2(x, 7, (j>k) was chosen arbitrarily, we have that d{y, fl{x, l3,ipk) < e 
for all a; € f2(x, 7, 0^). The same arguments show that S can be chosen small enough so 
that d{y,fl{x,j,4>k) < e for all x G i7(x, /?, -0^). Therefore it follows that 
dff(17(a;,/3,Vfc),17(x,7,0fc)) < e. □ 

Proposition 6.2. Lei /3 = {ipi, iPn}, /3' = {ip[, ip']^}, 7 = {0i, . . . , 0Ar}, anrf 
7' = {0']^, . . . , 0^} be orthonormal bases and suppose x G P7Y. T/ien /or any e > there 
exists a d > such that 

|Vol(f](x, A 7^,) n f)(a;, 7, 0fc)) - Vol(f](x, /?', V;-) n r!(x, 7', 0',))| < e 

for j,k = 1,2, ... ,N whenever — Vjll < ^ o.'f^d ||0j — 0^-|| < (5 /or = 1, 2, . . . , iV. 

Proof. Fix 77 > 0. Given an e > 0, there exists, by proposition l6.1l a (5 > such that HV'j — 
■0j|| < 5 and ||0j-0^|| < (5 for j = 1,2, . . ., TV implies that dH{Vl{x, (3,%l)j),Vl{x, (3',^'^)) < 
77 and d//(f^(x,7,0j),rj(a;,7',0^-)) < 77 . So 

n{x,[i,i^,) C B(r!(a:,/3',i/^;.),e) = r!(x, /?', t/^;.) U 
f)(a:,7,0fc) C B(f7(a:,7',0'fe),e) =l^(x,7',0fe)ur2 

where 

■.^B{n{x,P'A-).e)\n{x,P',i^'j) and B(r!(a:, 7', 0i.), e) \ f^(x, 7', 0fc) • 

Therefore 

n{x, /?, V,) n 17(2:, 7, 0fc) C n(a;, #■) n f)(a;, 7', 0',) U Tf U U 
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where := n{x, /?', V^ ) n T^, := n{x, 7', (f>'^) D T^, and := T^. Thus 

5 

It is clear that Vol(T/) ^ as e — > which shows that by choosing S small enough we 
get 

Voi(r!(.T, p, v^,) n Voi(r!(a;, 7, 0fe)) < /?', 7/.;) n n{x, 7', 0^)) + r? . 

Similar arguments show that by choosing S small enough we also have 

Voi{n{x, p', ^ij'j) n Voi{n{x, 7', (t>'k)) < ^i^, n{x, j,<Pk)) + v- 

Therefore we have established that 

|Vol(0(a;,/3',^;.) nO(ar,7',0'fe)) - Vol(0(ar,/3,^j) n^a;,7,<^fe))| <V- 

□ 

As before, let Sa := { ^4 e jCCH) \ = A} he the set of self-adjoint operators and 
V c Sa be the set of all self-adjoint operators with distinct eigenvalues. Note that V is 

open in Sa. 

Lemma 6.3. For each A £V there exist a neighborhood Na G V and smooth maps 

Uj : ^ K and 4>j : Ma ^ S{H) j = l,2,...,N 

such that for all B e Ma 

B<Pj {B) = aj {B)cf>j (B) j = l,2,...,N. 

Proof. Let ai < 02 < • • • < ajv be the distinct eigenvalues for A and Vi, • • • , V'JV the 
corresponding unit eigenvectors. Define a map 

by 

F^iB,^) := (B-a,(V,#)I)V'. 
This map is smooth and also satisfies F^{A, tpj) = O.Now 

B2F^{B, ip) ■ dilj = {B - aj{ipj\ilj)l )5ip - aj{ipj\6ip)ilj . 

Without loss of generality we can assume that Oj ^ because if it were then we could 
replace A be a nearby operator in V with aj 7^ 0. Expanding Sip as 

JV 

we get 

N 

B2F\B, ipj) ■ Sil) = ^ Cfe(afe - aj)'ij}k - ajCjipj ■ 
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Since aj ^ and ak—aj ^ for fc ^ j it is clear from the above result that {B, t{}j)T-i - 
Ti is an isomorphism. Therefore by the implicit function theorem there exists a neigh- 
borhood A/a of A and a map (l>j : Ma H such that 

F^{B, ^j{B)) = for aU B e A/'a and ^j{A) = . (6.11) 

Since WipjW = 1 we can by shrinking Ma if necessary define 



0, : ATa ^ S{n) : B 



\MB)\\ 
It follows from (|6.11|l that 

[B - a, (B) I (B) = for aU B E Ma 

where 

a,iB) -.^ a,{cj,,\4>,iB)) . 

This shows that aj : Ma ^ C is smooth. However, since B is self-adjoint, it has only 
real eigenvalues and we must have aj{MA) C M. □ 

The above lemma allows us to compare the joint distributions for the self-adjoint 
operators A,B and A',B' provided that A and A', and B and B' are close enough. 
To see this fix G V. Let Ma, Mb C P, (f>f : Ma S{H), (t>f : Mb ^ S{n), 
af : Ma ^ S{H), and af : Mb ^ S{H) be as in lemma lO Suppose A' e Ma and 
B' e A/fi, then we can define 

^f:^^f{A'), ^f:=0f(i?'), ijf:=^f{A'), ^f:=<^f(i?'), (6.12) 
a; af{A') , fo^ := a^(B') , a, a/(A) , 6, a^(B) . (6.13) 

For any x € P7i, we can compare the joint distributions Prob(y4 = aj, B = bk \ x) and 
Prob(yl' — a'j, B' = b'l^l x) by computing 

||Prob(A, B\x)- Prob(A', B' \ x)\\ 

N 

:= I Prob(A ^aj,B = bk\x)~ Prob(A' = a'j,B' = b',,\x)\. (6.14) 
j,fe=i 

Theorem 6.4. For any x E PH and e > there exists a S > such that if A, A' , B, B' E 

V and \\A- A'W, \\B - B'\\ < S then 

||Prob(A, B\x)- Prob(A', B' \ x)\\ < e . 

Proof. We know from the discussion following lemma 1^31 that for any 77 > there exists 
a 5 > such that if \\A - A'\\ < 6 and \\B - B'\\ < S then ||^/ - tJjflUf - V'f II < V 
for j — 1,2, . . . , N where ip^,^^ , -0^, and are as defined by (|6.12|l . Now, 

I Prob(A = cj, B = bk I x) - Prob(yl' = a'j,B' = b'^ \ x) \ 

= |Vol(r!(a:,/3-^,V'/)nl](a;,/3^,Vf)) - Vol(f](a:,/?^',^f )nr!(x,/5^',Vf' )) |. 
Proposition 16 . 21 then shows that if we choose 77 small enough then 

I Yo\{^i{x, /3^, ^f) n ^{x, /3^, )) - Voi(r!(x, p^' , ijf) n n{x, f}^' , v^f )) I < ^ 

and the proof follows immediately from (|6.14|l . □ 
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7 Conclusion 



We have presented a method for constructing joint distributions for an arbitrary number 
of observablcs on finite dimensional Hilbcrt spaces. In the case; where; the obscrvablcs 
commute, tlie joint distributions agree with the standard quantum meclianical ones. 
Furthermore, we have shown that that our joint distributions depend continuously on 
the observables and transform correctly under unitary transformations. 
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